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ABSTRACT 

We present upper limits on the X-ray emission for three neutron stars. For PSR J1840— 1419, with a 
characteristic age of 16.5 Myr, we calculate a blackbody temperature upper limit (at 99% confidence) 
of kT{£ < 24±\l eV, making this one of the coolest neutron stars known. PSRs J1814— 1744 and 
J1847— 0130 are both high magnetic field pulsars, with inferred surface dipole magnetic field strengths 
of 5.5 x 10 13 and 9.4 x 10 13 G, respectively. Our temperature upper limits for these stars are kT^ < 
123^33 eV and fcTJ°£ < 115^3 eV, showing that these high magnetic field pulsars are not significantly 
hotter than those with lower magnetic fields. Finally, we put these limits into context by summarizing 
all temperature measurements and limits for rotation-driven neutron stars. 
Subject headings: pulsar: general — stars: magnetic fields — X-rays: stars 



1. INTRODUCTION 



Measurements of the temperatures of neutron stars 
(NSs) are few and far between. They are often impos- 
sible as, for ages > 1 Myr, the star will have experi- 
enced significa nt cooling since its formation ( Yakovlev & 



Pethick 2004). The thermal emission from JNSs peaks 
in soft X-rays, below ~ 1 keV, and the luminosities are 
such that sources even a few kiloparsecs away are diffi- 
cult to detect. Currently, a total of 37 measurements, 



temperatures have bee 


n published (|Li et al. 2005 


Car- 


aveo et al.||2010| |Abdo 


2010; Spcagle et al. 2011 


c 


hang 


et ai. 


2011 


Zhu et al. 


2011 IMarelli et al. 2011 


Guvcr 


et al. 

'mm 1 


2012; 


Pires et al. 


2012 


). Of these, only 11 are for 



ing 5 millisecond pulsars for whom we have no meaning- 
ful estimate of age) , and only three of those greater than 
10 Myr: PSRs B0950+08, J2144-3933 and J0108-1431. 
The emission from PSR B0950+08, a 17.5 Myr pulsar, is 
non-thermal, with a power-law spectrum, and an upper 
limit of kTP? < 41 eV, or T < .48 MK, on the thermal 



contribution ( Becker et al.|20 04). To put this in context 
it is at the level ot the lowest meas ured temperature fo r 
any NS: kT{£ = 43 eV for Geminga ( |De Luca et al.|2005[ ). 
PSR J2144— 3933 is the slowest spinning pulsar known 
(P = 8.5 s) and an order of magnitude older (272 Myr) 
than PSR B0950+08 so that it is perhaps not surpris- 
ing that no X-ray emission has be en seen from it, with 



an upper limit of k T^ < 20 eV (Tiengo et al. 2011). 



Posselt et al. (2012 1 have recently reported the detection 
oi PSR JU10S=T531 (170 Myr) with fcT b °g = 110 eV, a 
source which seems to be undergoing a heating process. 

In a recent re-analysis of the Parkes Multi -beam Pulsar 
Surve y, PSR J1840— 1419 was discovered ( Keane et al. 
20101. This 6.6-s pulsar shows sporadic radio emission 
with strong single pulses detected, on average, every 
~ 10 rotation periods, at flux densities up to 1.7 Jy. 
At 16.5 Myr, PSR J1840-1419 has a very similar char- 
acteristic age to PSR B0950+08, and therefore might 



be expected to have somewhat similar X-ray proper- 
ties. This, plus its relative proximity at 850 pc, make 
PSR J1840— 1419 an excellent target for adding to our 
knowledge of the X-ray emission and temperature of old 
NSs. In this paper we present the result of a Chandra 
observation of PSR J1840-1419. 

One factor which has been suggested to be a major 
factor in NS cooling is the magnetic field strength. Sev- 
eral authors have suggested that NSs with higher mag- 
netic field strengths are hotte r than those with lower field 
strengths, f or the same age (Gonzalez et al. 2007 Zhu 

SHclT 



et al ||2009| |Qlausen et al.||20l0[ [Zhu et al.||20T^ _ 

an effect is theoretically predicteal Yakovlev fe Pethick 



2004[) and th e work of |Geppert et al.| ( |2004[ ) and |Perez 
Azorin et al. (2006) shows that high magnetic fields su 
press heat conductivity perpendicular to the field lines, 
naturally producing an anisotropic temperature distri- 
bution on the stellar surface with small hot regions at 
the magnetic poles. We examine this claim with XMM- 
Newton observations of two pulsars: J1814— 1744 and 
J 1847— 0130, both of which have high magnetic field 
strengths, of 5.5 x 10 13 and 9.4 x 10 13 G, respectively. 

Finally, we present an updated plot of fcXj°g vs charac- 
teristic age for all X-ray measurements of, and upper lim- 
its on, rotation-powered NS temperature made thus far, 
and re-examine the temperature-magnetic field strength 
relationship. 

2. OBSERVATIONS & RESULTS 

2.1. Calculating Limits 

For a non-detection, one could deter mine a count rate 
limit by following the procedure of Pivovaroff et al. 
(2000): (i) define an aperture on the image about the 
position of the source, and add up the counts in this 
region, N; (ii) define a background region of the same 
area (or appropriately scale to the same area) elsewhere 
on the image and add up the counts in the background 
region, B. The signal from the source is N — B and 
the noise is y/N + B. The signal-to-noise ratio is then 
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(N - B)/\/N + B. To set a '3 - a limit' count rate we 
can then solve for N in 9(N + B) = (N - B) 2 , where 
B is known. We can see that in the A ^> B 'photon- 
rich' case this tends to the expected Poisson value in 
the absence of background sky noise. However, in the 
'photon-poor' case we clearly can not use this expres- 
sion, e.g. for B — 0.35, N = 10 naively suggests a 99.7% 
detection/limit, however we know that the Poisson prob- 
ability for 10 counts due to noise is ~ 10 _13 %; for such a 
background rate 3 counts represents a 99.8% limit. The 
observations presented here are in the photon -poor sce- 
nario, so that using the Pivovaroff et al. ( 2000 ) procedure 
would result in a much poorer limit than what has truly 
been obtained. 

To convert count rate limits to flux limits we need ei- 
ther some kind of absolute scale, e.g. from observing 
a standard source of known flux, or a model for our 



sourc e, e.g. an assumed spectral form (McLaughlin et al 
2003 1 . Below we calculate one set of limits for a black- 
body model (whence we obtain temperature limits as a 
function of the emitting radius), and another set of limits 
for a non-thermal model with a power law with negative 
index of 2.0. 

2.2. PSR J1840-U19 

PSR J1840-1419 is an old pulsar (16.5 Myr) with a 
6.6-s spin period. Due to its sporadic radio emission it 
was detected in a s earch for single pu lses, rather than in a 
periodicity search ( Keane et al.|2010 ). Pulsars discovered 
in this way are of ten referred to as "RRATs" ( |Keane fc 
McLaughlin|2011[ ). The properties of PSR J1840-1419, 
both measured and derived, are given in Table [l] 

On 2011-02-20, we performed a 10-ks observation of 
PSR J1840-1419, using the ACIS-S detector on Chandra 
in the energy range 0.2 — 10.0 keV. We detected only one 
photon within the ~ 1 arcsec (3 pixe l) error circle of th e 
position derived from pulsar timing ( Keane et al.||20lT ). 
Assuming that this 1 count is consistent with background 
noise (the background rate being 0.24 counts/pixel), and 
ignoring the Poisson nature of the source, the count rate 
limit is 3 x 10~ 4 counts/s (at 99% confidence). 

In the case of a blackbody model this count 
rate limit implies a temperature limit of kT££ < 

where the dependence on 



-0.39 



24^eV(i?£° b /10 km) i 
i?^ is a purely empirical fit to a simple power law of 
the form T\qR" , where Rio is the blackbody radius 
in units of 10 km and T10 is the temperature when 
i?io = 1: see top left panel of Figure [I] The error 
bars reflect both NE2001 ' worst case' errors of a fac- 
tor of two in the distance ( |Cordes fc Lazio] 2002 1 as 
well as factor of two uncertainties in the neutral hy- 
drgoen column density (see below). The bottom right 
panel of Figure [l] shows all measured neutron star tem- 
peratures as a function of age. The PSR J1840— 1419 
limit is amongst the coolest of all published limits. 
The corresponding flux and luminosity limits are /bb < 



5.0 x 10- 14 (i?^/10 km) 



11 erg s 1 cm 2 



and < 

4.3 x 10 30 (i?j^/10 km) 161 erg s"\ respectively. In the 
case of the non-thermal model (power law, with negative 
index of 2.0), the flux and luminosity limits are f nt < 
3.0 x 10~ 15 erg s" 1 cm" 



2 and < 2.6 x 10 29 erg s" 



respectively. 

The above estimates use a neutral hydrogen column 



density of Ah = 6 x 10 20 cm~ 2 which is derived from 
the dispersion measure, assuming 10 neutral hy drogen 
atoms per free electron (Seon & Edelstein 1998). Pre- 
dicted count levels were calculated using the standard 
PIMMS (Portable, Interactive Multi-Mission Simulator) 
tooQ 

2.3. PSRs J1814-1744 & J1847-0130 

We also present the results of two XMM-Newton 
observations using the PN detector and medium filter in 
the energy range 0.15 - 15.0 keV. PSR J1814-1744 was 
observed for 6.1 ks on 2004-10-21, and PSR J1847-0130 
was observed for 17.0 ks on 2004-09-14. The observa- 
tions of these young, high-B pulsars also resulted in 
non-detections. As above for PSR J1840— 1419 we derive 
temperature upper limits for these sources of: kT^ < 

for PSR J1814-0130, 



123t|°eV(i2ft/10 km)- 



0.23 



and kT t 



bb 



< 



115 



+ 16 
-33' 



:V(i?g° b /10 km)" 



for 



PSR J1847-0130 (see Figure [TJ. These lim- 
its are much less constraining given the much 
larger estimated distances (7.7 and 9.8 kpc re- 
spectively). The associated flux and luminos- 
ity limits for the blackbody scenario are: /bb < 



\1.77 



2.6 x 10- 13 (i?g|/10 km) 

1^/10 km) 177 erg s 



erg s 1 cm 2 and Lgj^ < 



3.0 x 10 33 (i?^/10 km) 177 erg s" 1 for PSR J1814-1744; 
and /bb < 3.2 x 10- 13 (.Rj£/10 km) 179 erg s" 1 cm~ 2 
and L^ h < 2.3 x 10 33 (i?^,/10 km) 179 erg s" 1 for 
PSR J1847— 0130. 

The non-thermal limits are: / nt < 8.9 x 
10~ 15 erg s _1 cm -2 and L££ < 1.0 x 10 32 erg s _1 for 
PSR J1814-1744; and: / nt < 5.1 x 10~ 15 erg s" 1 cm~ 2 
and L™ < 3.6 x 10 31 erg s" 1 for PSR J1847-0130. 
For PSR J1814-1744 the value of N H derived from the 
dispersion measure was ~ 50% higher than the maxi- 
mum Galactic value for this line of sight, as calculated 
by the standard COLDEN tooQ so the COLDEN value 
of Nn = 1.6 x 10 22 cm~ 2 was used. For this reason ATh 
is unlikely to be underestimated, but we have allowed 
for the possibility that it may be overestimated by as 
much as a factor of two. Likewise, for PSR J1847— 0130, 
the neutral hydrogen column density was derived to be 
Ah = 2.1 x 10 22 cm~ 2 . This value of Ah is very close to 
the maximum value for this line of sight, so that it too 
is unlikely to be an overestimate. 

3. DISCUSSION & CONCLUSIONS 

Although this paper reports three null results we deem 
it to be of the utmost importance to report such inves- 
tigations to avoid duplicated efforts, wasted telescope 
resources, etc. Table [T] summarises the temperature, 
flux and luminosity limits for the sources reported here. 
From the non-thermal luminosity limits we can deter- 
mine upper-limits on the non-thermal X-ray efficiency, 
V = L nt /E, where E is the spin-down energy loss rate. 
Although the luminosity limits for the XMM-Newton ob- 
servations are ~ 2 orders of magnitude poorer than for 
the Chandra observation of J1840— 1419, the E values 
are correspondingly higher, s uch that for all 3 puls ars the 
limit turns out to be 77 < 0.2. Possenti et al. ( 2002 ) deter- 
mined, for their study of 39 pulsars, that all sources had a 

1 http:/ /asc. harvard.edu/toolkit/pimms.jsp 

2 http:/ /asc. harvard. cdu/toolkit/colden.jsp 
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value of X] less than ry max = 10 _18 ' 5 (£'/erg s -1 ) 48 . Our 
r] limits are all significantly higher than this predicted 
value. However, we note this critical 77 value is arrived at 
upon considering pulsars studied in the 2 — 10 keV energy 
range, narrower than the observations reported here. 

These observations bring the total number of rotation- 
powered NSs with detections of thermal emission, or up- 
per limits thereupon, to 49. In Figure[l]we show our lim- 
its, along wi th all previous measurements and limits. |Zhu 
et al. (2011) suggest that there is a "hint" that the high 
magnetic field pulsars are hotter than the low magnetic 
field pulsars. To quantify this we first ignored the up- 
per limits and divided the remaining data into 'low' and 
'high' magnetic-field-strength sources, below and above 
an arbitrarily chosen field strength value of 10 13 G. We 
also excluded three 'low' magnetic-field-strength sources 
as their fitted blackbody radii are so much smaller than 
those of the rest of the sources (at 33, 43 and 120 m) so 
that their therm al emission is thought t o be due to some 
heating process ( Misanovic et al.||2008 ). Comparing the 
two distributions with a Kolmogorov-Smirnov (K-S) test 
we find a K-S statistic of D — 0.33, and thus a proba- 
bility of 0.31 that these two distributions are the same. 
There is therefore no evidence that the high magnetic 
field pulsars are hotter than those with lower magnetic 
fields. However we note the caveat that, in addition to 
the age estimates being uncertain (as mentioned above) , 
the magnetic field estimates are uncertain, and should 
only be considered accurate to within an order of mag- 
nitude. Knowledge of the inclination angles between the 
magnetic and rotation axes would be needed to estimate 
the magnet i c field strengths more accurately (e.g. as in 
Spitkovsky| (|2006[)). 
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Table 1 



Measured and derived properties, and limits for 3 pulsars. 


Quantity 


J1840— 1419 


J 1814— 1744 


J1847— 0130 


P (s) 


6.6 


4.0 


6.7 


P (is Is) 


6.3 


745 


1270 


T c (kyr) 


16500 


85 


83 


B (10 12 G) 


6.5 


55.1 


93.6 


E (10 30 erg s" 1 ) 


l.O 


468 


167 


DM (cm" 3 pc) 


19 


792 


667 


D (kpc) 


0.9 


9.8 


7.7 


Telescope 


Chandra 


XMM 


XMM 


Instrument 


ACIS-S 


PN 


PN 


Date 


201 1 -02-20 


2004-1 0-21 


2004-09-1 4 

£* UUi VJZ7 -Lt: 


Energy Range (keV) 


0.2 - 10 


0.15 - 15 


0.15 - 15 


T obs (ks) 


10.0 


6.1 


17.0 


Nh (10 20 cm" 2 ) 


6 


157 


205 


kT™ (eV)f 


<U±\l 


< 123 +2 ° 


< 115+ 33 


T°° (MK)t 


< 0.28«;i 9 2 


< 1 42+ - 22 




/bb (IO" 14 erg cm" 2 s" 1 )! 


< 5.0 


< 6.1 


< 8.9 


L bb (10 30 ergs- 1 )! 


< 4.3 


< 692 


< 640 


/at (10 15 erg cm 2 s x ) 


< 3.0 


< 8.9 


< 5.1 


L nt (10 30 ergs" 1 ) 


< 0.26 


< 102 


< 36 


7? = L n t/S 


< 0.26 


< 0.22 


< 0.22 



Note. — The measured and derived properties for the three pulsars discussed in this paper. The spin period, P, and its derivative, 
P, are derived from pulsar timing techniques, and r, B and E, the characteristic age, magnetic field strength and spin-down energy 
respectively, are inferred from these jK eane et a l. 2011). The distance, D, is derived from the measured DM using the NE2001 model for 
the Galactic free electron distribution (Cordes & Lazio 2002). The 99%-confidence temperature limits are derived as outlined in the text. 
Non-thermal flux limits are given for a power law with negative index of 2.0. f these limits are for an emitting radius of 10 km, and scale 
as (R??/10 km)°, where a equals —0.39, —0.23 and —0.21, respectively for the three sources (see main text). \ The blackbody flux and 
luminosity limits similarly scale as (R^/10 km) 2-0 . 
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10 10 10 10 10 

Characteristic Age (yr) 

Figure 1. The first three plots show the temperature limit (i.e. which produces the observed count rate limit for an assumed 
blackbody model) as a function of the assumed emitting radius, for the Chandra observation of PSR J1840— 1419 (10 ks) and 
the XMM-Newton observations of PSRs J1814-1744 (6 ks) and J1847-0130 (17 ks). The points are fit with a simple power-law 
of the form Tio(_Rio) a . In each case the middle curve is for the n ominal dispersion-mea sure-derived distance. The cooler and 
hotter curves are for 'worst-case' distance errors of a factor of 2 ( Cordes fc L azio 20021. For clarity of presentation we show 



here only the curves for the nominal Nh values quoted in Table [T] The fourth plot shows all known temperature values as a 
function of characteristic age. Sources with inferred magnetic field strengths > 10 13 G are marked in blue, whereas those with 



lower values are marked in black. Open circles denote upper limits, 
with large circles for clarity. Note that the symbols for PSRs J1814- 
and derived limits are very similar. 



The three upper limit presented in this paper are marked 
1744 and J1847— 0130 are largely overlapping as their ages 



